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ABSTRACT 

The turbulence equations are closed by s|x»cification of initial condi- 
tions (usinfc, either a Taylor or an exponential series) and by a modified 
Kovasz nay -type closure Good results for large times are obtained only 
for the initial-conditions closure used with four or more terms of an 
exponential series. The evolution of all of the initially-specified spectra 
can lx* calculated rather well from the theory. From a fundamental 
standpoint the method thus seems to lx' satisfactory. 
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I. INTRODUCTION 

From a practical standpoint it would be advantageous to be able to 
calculate the evolution of the turbulent energy spectrum by specifying 
only that quantity initially. Unfortunately, because of the coupling be- 
tween the members of the infinite hierarchy of multiixiint correlation or 
spectral equations, it apix'ars that we are not able to do that, so that a 
satisfactory theory would seem to require the initial specification of a 
number of interacting quantities. The prediction of the evolution of the 
energy spectrum in fact requires the specification of an infinite number 
of initial multipoint correlations or spectra (or functionals of these quan- 
tities). * Those quantities, together with the correlation or spectral equa- 
tions, can be used to calculate the initial time derivatives of the correla- 
tions or spectra. The evolution of the turbulence can then be obtained by 
using an exponential time series which is an iterative solution of the Navier- 
Stokes equations. When the problem is formulated in this way, the correla- 
tion or spectral equations are closed by the initial specification of the tur- 
bulence and no assumption is necessary for closing those equations. 

Of course in practice we can specify only a finite number of the lower- 
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order quantities. This has been called the "gap problem." It is the 
problem of bridging the gap between the infinite number of correlations 
which would theoretically be necessary for calculating the evolution of 
the turbulence, and a finite specifiable number of correlations. Most 
workers have attempted to bridge the gap by assuming that the initial 
distribution of turbulent fluctuations is exactly Gaussian (zero odd-order 
correlations). However, that is an artificial initial condition, probably 
never realized for real turbulence. The importance of accurate initial 
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conditions is shown, for instance, by the data of Ling and Saad. where 
measurements were made downstream from a turbulence-producing 
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waterfall. The turbulence decay law for the Initial conditions produced 
by the waterfall is considerably different from that for initial conditions 
produced by a grid. As will be seen later, the skewness factors for the 
velocity gradient are also different. 

Here, as in Ref. 1, we bridge the gap simply by specifying a suffi- 
cient number of initial correlations or their spectral equivalents to cal- 
culate the evolution accurately. Fortunately, we do not ha\e to specify 
the multipoint correlations or their spectral equivalents themselves 
(those would lx* extremely difficult to measure) but only two -point func- 
tionals of the multipoint spectral quantities. ^ It will be seen that the 
evolution of all the quantities which are specified initially can be cal- 
culated. If, on tiie other hand, a large number of initial conditions were 
specified in order to predict the evolution of say one quantity, it might 
be objected that the initial conditions were chosen to make the theory 
agree with experiment for that one quantity. However, that objection 
cannot l>e made if the evolution of all of the quantities which are specified 
initially can be calculated, as in the present theory. From a fundamental 
standpoint the calculation of the evolution of those quantities is all that 
might be expected from a theory of evolving turbulence. 

II. ANALYSIS 

In Refs. 1 and 4 we gave the basic theory for closure by specification 
of initial conditions and calculated results for low and moderate turbulence 
Reynolds numbers (1I X between 3 and 70, where R^ is the Reynolds num- 
ber based on the Taylor microscale and the root-mean-square turbulent 
velocity fluctuation). Here we compare calculated results with the higher 
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Reynolds number data of Ling and Saad. The Reynolds numbers for those 
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data were high enough to obtain a -5/3-power region in the energy sjjec- 
trum (H x between 300 and HOO). 

The exponential-series expression for the energy spectrum function 
E was obtained in Ref. 1 as 

EU.t) -- B 2 U)exp[-2i'K 2 (t - tj)] + ^ ' B 2 (K)exp[-2b § (x)(t - t,)], (1) 
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where k is the wave number, t is the time, tj is the intial time, and v 
is the kinematic viscosity Equation (1) gives the evolution in time of the 
energy spectrum from an initial state which is specified by the B’s and 
b’s in the equation. The first term is the well-known expression for the 
decay of E in the final period (weak turbulence approximation). The 
rest of the terms in Eq. (1) therefore give the contribution of inertia to E. 
In the present note we will retain a maximum of four exponential terms 
in Eq. (1). This is one more term than it was necessary to retain for the 
low and moderate Reynolds number data considered in Refs. 1 and 4. 

With four terms retained, we will have to specify seven spectra at tj in 
order to evaluate the functions B, Bj, B 0 , Bg, and bj, b 0 , and bg. 
Evidently we need more spectra to describe the initial turbulence at higher 
Reynolds numbers because a wider range of eddy sizes is excited, and the 
turbulence structure is more complicated than at lower Reynolds numbers. 
The spectrum functions, in addition to E(x, t), which we will specify ini- 
tially are T(k.I), VU,t), R(k,0, S(k , t) . L(k. t), and M(k, 0, which 
are given by 


T = 0E/0t + 2 pk 2 E, 

(2) 

V = ?T/0t + 2 ck 2 T, 

(3) 

R = av/?t + 2 i*k 2 V, 

(4) 
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S = PR/ at f 2t>K 2 R, 

(5) 

L 38/ at + 2ck 2 S, 

(6) 

M = PL/ Pt + 2 ck 2 L. 

(7) 


The quantity T is the well-known energy -transfer function, and V, R, 

S, L, and M are, respectively, two-point functionals of three-, four-, 
five-, six-, and seven-point spectral quantities. ' The latter five quan- 
tities are somewhat similar to T, inasmuch as Uiey contain the effects 
of transfer between wave numbers of T, V, R, S, or L. However they 
differ, in that they also contain other effects, so that the areas under 
those spectra are not necessarily zero, as in the case of T. Equations (2) 
through (7) are obtained, respectively, from two-, three-, and four-, five-, 
six-, and seven -point spectral equations. 1 They c; n be thought of as two- 
point alternatives to the multipoint spectral equations and are much easier 
to work with. 
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For comparing the theory with the experiment of Ling and Saad, the 
experimental input can be conveniently obtained from an empirical equa- 
tion for E, Eq. (8) in their paper. The higher-order spectra were not 
measured directly in their experiment, but could be calculated from their 
equation for E and Eqs. (2) through (7). Except for experimental error 
those values will be the same as those that might have been measured 
directly. The B’s and b’s in Eq. (1) can be related to the initial spec- 
tra by successive differentiations of that equation with respect to time, 
setting t tj, and using Eqs. (2) through (7). 

The use of Ling and Saad’s Eq. (8) for obtaining the initial conditions 
might be taken as ar. indication that it is necessary to know the decay 
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law to predict it. However that is not the case. Their equation was 
used because the initial spectra could not be measured directly in their 
experiment. Moreover we do not need to know E for the whole decay 
period, but only at enough early times to calculate the required initial 
spectra, the latter being a much smaller amount of information. At anv 
rate the nature of the problem seems to be such that we do not have a 
reasonable alternative to the procedure followed here. 

III. RESULTS AND DISCUSSION 
Before giving results obtained from Eq. (1), we will consider a 
Taylor series with a maximum of seven initial spectra (the same maxi- 
mum number that will be used with Eq. (1)), and a modified Kovasznay- 
type closure (modified to include an effect of initial T).^ Results for 
those calculations are given in Fig. 1, where the initial state is speci- 
fied at tj = 0.0048. The quantity u. is a velocity component, the over- 
bar indicates an averaged value, and the stars indicate that quantities 
have been nondimensionalized by the kinematic viscosity and an experi- 
mental constant A (the proportionality constant in the power decay law 
for Uj) having the dimensions ^(length) 1 * (time) 1 The turbulent 
energy ( 1 / 2 ) 11 ^ is obtained by integrating E over all wave numbers. 
Figure 1 indicates that the Taylor-series results agree with experiment 
only for small times. Evidently, many more terms (and initial spectra) 
would be required in order to obtain accurate results for u^ at large 
times. Note that even if enough terms were x-etained in the Taylor series 
to give accurate results for u^, the decay of the higher-order spectra 
which would then have to be specified initially could not be accurately 
calculated. Thus the use of a Taylor series in the present problem does 
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not give a satisfactory solution, regardless of the number of terms re- 
tained. 

The modified Kovusznay -type closure is in somewhat better agree- 
ment with experiment than the Taylor series, but at large times the agree 

ment is still not good. This is in contrast to its good agreement at moder 
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ate and low Reynolds numbers. It was introduced in Ref. 5 in an effort 
to reduce the required number of initial spectra. Evidently that effort is 
successful only for moderate and low Reynolds numbers. It is possible, 
of course, that a more sophisticated method might be more successful 
(e. g. see Ref. 6). 

A comparison between theory and experiment using the exponential 
series (Eq. (1)) is given in Figs. 2 to 8, where the initial state is again 
specified at tj = 0.0048. As in Refs. 1 and 4. unphysical singularities 
occasionally occurred in the theoretical spectra. Inasmuch as the un- 
physical values were localized, particularly in the higher-order approxi- 
mations, they were simply omitted in plotting the curves. 

Figure 2 gives a comparison between theory and experiment for the 
decay of turbulent energy. Theoretical curves are shown for one, two, 
three, and four exponential terms retained in Eq. (1). The curve for 
four terms is in good agreement with experiment for the whole decay 
period. Comparison of the curve in Fig. 1 for one term retained (weak 
turbulence approximation) with the experimental curve shows the effect 
of inertia on the decay process. In contrast to the results for moderate 
Reynolds number, 4 where inertia and viscous effects were of the 
same order of magnitude, the inertia effects for the present higher 
Reynolds number results are at least an order of magnitude greater than 
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the viscous effects. Thus most of the decay at high Reynolds numbers 
is due to inertial self -inter twtion of the turbulence, rather than to viscous 
effects. Figures 3 and 4 shov how the energy and the transfer spectra 
decay with time. 

Figures 5 and 6 compare theory and experiment at a late time for 
all of the spectra which are initially specified to describe the initial tur- 
bulence. The prediction of the decay of all seven of the spectra which 
are specified initially is rather good. Thus the present theory appears 
to solve an initial-value problem for higher Reynolds number turbulence 

in which the decay of seven initially-specified spectra is predicted. 
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Although the initial dissipation spectrum k E is not specified in- 
dependently, because of its importance in turbulence theory it is com- 
pared with experiment and with the energy spectrum at a late time in 
Fig. 7. Again, good agreement is indicated. The separation of the 
energy and dissipation spectra is good evidence that we are dealing with 
high Reynolds number turbulence. 

Another important quantity is the skewness factor of the velocity 
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gradient. That quantity can be calculated from the spectra of E and T. 
The plot in Fig. 8 indicates excellent agreement between theory and ex- 
periment for the skewness factor. The difference between the trend in 
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Fig. 8 and that in grid-generated turbulence is probably due to the differ- 
ence between initial conditions for grid-generated turbulence and the 
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waterfall-generated turbulence considered here. 
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Figure 7. ■ Comparison of theoretical energy dissipation 
spectrum it l lit* time it* • 0.01 1 with experiment of 
reference 3 and with energy spectre. 
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Figure 8. - Comparison of evolution of theo- 
retical skewness factor of velocity gradient 
with experiment of reference 3. 



